Journal of Medical Entomology

Larval source management (LSM) is a complementary malaria control intervention, yet evidence to guide context-specific implementation remains limited. Nigeria's recent national commitment to LSM scale-up makes the need for operational evidence particularly urgent. Informal settlements embedded within wards of differing dominant settlement archetypes may present distinct Anopheles larval habitat profiles with implications for how LSM strategies should be tailored. We evaluated Anopheles larval habitats within informal settlement areas across wards with contrasting settlement archetypes in Ibadan metropolis, Nigeria, to inform targeted larval source management. Potential breeding habitats were surveyed in dry and wet seasons within informal settlement areas across three wards -- Olopomewa, Challenge, and Agugu -- representing formal, informal, and slum settlement-dominant archetypes respectively. Habitats were characterized and assessed for Anopheles larval presence. Pareto analysis identified habitats accounting for 80% of larval abundance. Breeding habitat density per km{superscript 2} was estimated using a simulated pathway technique. Associations between mosquito dispersal scale and household malaria infections identified through Rapid Diagnostic Testing were evaluated using kernel-based distance-decay weighting. Environmental drivers of habitat suitability were modeled in MaxEnt. Of 420 potential breeding habitats identified, 31 (7.4%) contained Anopheles larvae, predominantly during the wet season (26, 83.9%). Puddles, dug wells, drainages/gutters/ditches and canals accounted for 80% of site-level larval abundance when standardized by sampling effort. Larval and breeding habitat density were highest in Agugu, the slum-dominant ward, across both seasons. Modeled mosquito dispersal scale showed best fit at 30-32m in Challenge (OR 1.41, 95% CI: 1.05-1.89) during the wet season and 16-18m in Agugu (OR 1.29, 95% CI: 1.04-1.60) during the dry season. Habitat suitability in Agugu was higher farther from large water bodies and in areas with higher population density and positive Normalized Difference Water Index values. In Challenge, suitability was higher in areas with lower nighttime light levels, positive Normalized Difference Water Index values, and negative Normalized Difference Moisture Index values. Further studies incorporating multiple wards across diverse urban settings are needed to determine whether differences in larval ecology between settlement archetypes provide a reliable basis for planning larval source management.

Background Dengue is intensifying globally due to climate change, urbanization, and land use changes. In Nepal, dengue has expanded from lowland regions to higher altitudes, with record outbreaks in 2022 and 2023. However, reliance on passive surveillance and hospital-based studies may underestimate community level infection burden. Methods We conducted a population-based serologic cohort study in Kathmandu and Kavrepalanchok districts, Nepal, enrolling a geographically representative, age stratified random sample of residents aged 0 to 25 years from pre-defined hospital catchment areas. Enrollment occurred in two phases: Phase I (February 2019 to April 2021) with follow-up visits at approximately 3, 6, and 12 months, and Phase II (February to June 2023) revisiting the original cohort. At each household visit, we collected capillary blood samples by finger-prick onto filter paper and tested the samples for IgG responses against dengue-derived recombinant antigen using InBios DENV DetectTM ELISA. Serostatus was classified using the manufacturer's recommended immune status ratio (ISR) cutoffs. We calculated seroprevalence at each time point and estimated seroincidence rates by identifying seroconversion events per 1,000 person-years. We assessed risk factors using multivariable regression models. Results Between 2019 and 2023, we enrolled 840 participants and collected 2,082 blood samples. The overall seroincidence rate was 33.8 per 1,000 person-years (95% CI [24.9 to 45.0]), with the highest rates in urban Kathmandu ([105.7], 95% CI [75.1 to 144.4]). Seroincidence increased with age and over time from 46.1 in 2019 to 51.0 in 2023. Participants living with a dengue-positive individual in the same household (adjusted RR [4.65], 95% CI [2.72 to 8.0]) and households with water-filled flower basins (adjusted RR [2.53], 95% CI [1.28 to 5.74]) had significantly higher risk of seroconversion. Conclusions This study reveals a significant and increasing burden of dengue infection in the Kathmandu Valley between 2019 and 2023. highlighting an urgent need for immediate public health interventions to mitigate dengue's rise in Nepal's higher-altitude regions.

The World Health Organization (WHO) validated Sri Lanka had eliminated lymphatic filariasis as a public health problem in 2016, the second country in Southeast Asia to attain this status. However, post-validation surveillance has identified sporadic cases of brugian filariasis. The reemergence of Brugia malayi infections in Sri Lanka warrants urgent investigations. Recent studies have shown that the parasite responsible for the reemergence is a novel zoonotic Brugia sp. maintained among dogs that is closely related but distinct to the human-infecting B. malayi species. The current study employed morphological and morphometric assessments, revealing that this novel zoonotic Brugia sp. is within the B. malayi morphological range. Molecular characterization of three genomic regions, the nuclear genomic region SLXI, the non-coding region HhaI, and the mitochondrial genomic region COXI confirmed it as a genetic variant more closely related to B. malayi than to B. pahangi. Phylogenetic analysis further indicated it as a distinct genomic variant, closely related to a B. malayi-like parasite reported from India. Notably, that same parasite was identified in infected humans, animals, and potential vector mosquitoes. This, together with the detection of both human and animal blood within the same brugian infective mosquitoes, and delineating the canine origin of the parasites in human infections, provides compelling evidence supporting zoonotic transmission of this parasite. To our knowledge, this is the first report demonstrating the presence of the same brugian parasite in humans, domestic animals, and potentially infective mosquitoes in Sri Lanka, supported by multi-genomic evidence. The recent identification of multiple potential mosquito vector species suggests that this parasite may have undergone adaptive changes, facilitating its ability to overcome the species barrier. These findings substantiate the long-held hypothesis of zoonotic transmission of the reemerged brugian parasite, highlighting significant implications for ongoing surveillance and control strategies.

Dengue is a mosquito-borne, viral disease of increasing public health significance. Currently, most public health interventions target the vector, with efficacy dependent on timing within the season. Whilst seasonal profiles have been characterised in some endemic settings a global assessment is lacking. Here, we develop and apply a proportion-based measure of dengue seasonality to reported case time series from 1990 to 2024 across 106 countries and territories, the largest assessment of this phenomenon to date. We identify regional differences in seasonality such that every month of the year saw cases peak in at least one country or territory. Latitude was identified as influencing seasonality, with cases peaking between March and April in the southern hemisphere and July and October in the northern hemisphere. Equatorial locations displayed flat seasonality, and amplitude increased with distance from the equator. K-means clustering identified three seasonal profile types: two with pronounced seasonal outbreaks (with distinct peak timing and shape) and one with flatter, more endemic transmission. Peak month timing covaried among locations within the same seasonality cluster, with phase differences meaning that information on shifts in peak timing may be available several months in advance in some settings, of potential significance for prediction and intervention planning. Beyond aiding public health planning, identification of seasonal clusters suggests that information on dynamics in one location could be leveraged to improve forecasting power in others with similar seasonal dynamics.

Malaria remains a severe health problem in endemic regions because people lack adequate diagnostic tools, leading to delayed medical care and elevated death rates. This research introduces a dual-mode artificial intelligence system that uses two complementary models to enhance malaria pre-screening and diagnosis. The patient-centered model uses multivariate logistic regression to analyze biosignals, including heart rate, body temperature, and oxygen saturation, collected through a wearable sensor prototype and a mobile interface for symptom analysis. The system enables patients to begin self-assessment to determine their level of need before scheduling a doctor's appointment. The clinician-centered model represents a customized convolutional neural network that uses annotated microscopy images of red blood cells to achieve 94.84% accuracy, 95.71% precision, 93.87% recall, 94.78% F1 score, and 0.84 Area Under Curve (AUC). The patient model achieved 94.6% accuracy and an AUC of 0.985 using a 70/30 train-test split. These systems work together to create a layered diagnostic system that can operate independently or together to detect malaria at an early stage, especially in areas with limited resources. The findings demonstrate that wearable biosignal data integration with image-based deep learning can produce dependable, scalable, and user-friendly systems for malaria pre-screening. Keywords - malaria diagnosis, artificial intelligence (AI), convolutional neural networks (CNN), wearable biosensors, multivariate logistic regression

Introduction: Uganda has the highest number of Malaria cases in East and Southern Africa and accounted for 3% of deaths in 2020. This calls for prevention, early diagnosis and treatment of Malaria most especially among children whose condition can progress to severe Malaria within 24 hours. While numerous interventions have been put into place to prevent malaria transmission, delays in diagnosis and treatment of Malaria when ill can lead to further mortality. Therefore, factors associated with delayed access of care among children under five with malaria and their outcomes need to be explored. Methods: A cross sectional study was carried out. The target population was parents/caretakers to children under five with malaria at Mbale Regional Referral Hospital. A consecutive sampling technique was used on the target population. Quantitative data was collected using researcher administered questionnaires designed consistent with the research objectives. Collected data was analyzed using STATA version 15. Results: Among the 216 children under five admitted at Mbale regional referral hospital with Malaria, 59.26% received care from a health facility 24 hours after symptom onset. The most significant predictors of delay in seeking care were the caregiver/ parent having attained tertiary education (AOR=7.1, p value=0.02) and initially implementing other measures other than giving medication/herbs before taking a child to the health center (AOR=4.1, p value=0.00). Conclusion: Despite the numerous interventions put into place to curb the spread of malaria and to manage malaria, delayed access of care remains a significant contributor to the adverse effects of malaria among children under five. Health education on the impact of delayed access of care should be intensified at all levels of healthcare.

Each year, over 50 million children receive preventive malaria treatment. However, to date there has been no consensus on the most effective antimalarial drugs to use, especially given geographic differences in drug resistance. Here, we conduct a systematic review comparing the effectiveness of the most commonly used antimalarial chemopreventive regimen, sulfadoxine-pyrimethamine plus amodiaquine (SP+AQ), with other antimalarial drugs in preventing new infections. We searched MEDLINE, Embase, Global Health, PubMed and WWARN clinical trial databases until 06 December 2025 for studies satisfying the inclusion criteria. Studies were included if they were peer-reviewed, randomised-controlled studies in Africa, measuring incidence of infection or clinical episodes of Plasmodium falciparum malaria for at least 28 days post-treatment. We also compiled data on the prevalence of markers of resistance in the parasite dhfr, dhps and mdr1 genes in the study areas. We conducted meta-analyses of incidence rates, with subgroup analyses by drug resistance levels. This review is registered on PROSPERO (CRD42024577149). We identified 27 studies representing 38,252 participants in 32 sites across 13 countries. In pooled analysis, SP+AQ reduced incidence of malaria by 54.6% (95% CI: 33.8-68.8%) compared to SP alone, including significantly outperforming SP even in areas with low SP resistance. These findings suggest that countries currently using SP alone for chemoprevention should consider switching to SP+AQ. Where AQ resistance remains low, available evidence suggests SP+AQ remains efficacious for malaria chemoprevention. SP+AQ was comparable to the artemisinin-based treatment, dihydroartemisinin-piperaquine across all studies (incidence rate ratio 0.93; 95% CI 0.78-1.11). By resistance levels, SP+AQ had slightly higher efficacy in areas with low SP and AQ resistance but had comparable or slightly lower efficacy in areas with higher resistance. Using artemisinin-based treatments for chemoprevention must be balanced against the risk of worsening artemisinin resistance in Eastern and Southern Africa. This study was funded by the UK Royal Society.

An epidemic's expected course is determined by the magnitude and timing of a typical person's infectiousness --- captured, in turn, by the basic reproduction number and the generation-time distribution. These fundamental, population-average quantities can mask individual-level variation that shapes how an epidemic actually unfolds: for example, individual variation in the magnitude of infectiousness (overdispersion) creates superspreading, a key feature of the SARS-CoV-1 and SARS-CoV-2 epidemics. However, the impact of individual variation in infectiousness timing is less well understood. Here, we demonstrate that individual infectiousness timing varies substantially and to different degrees across pathogens. For some common pathogens, including influenza, measles, and SARS-CoV-2, infectiousness is "bursty", or highly concentrated and variably-timed across individuals: for example, the window of appreciable infectiousness for SARS-CoV-2 may last for roughly a day, vs. the 9--12 days usually quoted. We show that bursty infectiousness creates superspreading without inherent superspreaders, makes epidemic timing more variable, amplifies the time-sensitivity of common interventions, and complicates inference of key epidemiological parameters. Together with the reproduction number, the generation-time distribution, and overdispersion, burstiness completes a family of basic parameters that govern how epidemics unfold.

Background: Opioid use disorder (OUD) is heritable, yet most genome-wide association studies (GWAS) have focused on European populations, leaving the genetic architecture of OUD in non-European populations underexplored. Methods: We conducted GWAS of OUD across three ancestries using electronic health records and genomic data from 52,357 All of Us Research Program participants (8,912 cases; 43,445 matched opioid-exposed controls; 48.5% female). Participants were stratified into European (EUR), African (AFR), and Admixed American (AMR) ancestry groups for logistic regression GWAS, with independent replication in the Million Veteran Program. We then applied the deep-learning model AlphaGenome to predict the tissue-specific transcriptomic and splicing consequences of top risk variants across 13 reward-pathway brain regions. Results: We identified and replicated a novel DDX6 risk locus, alongside established OPRM1 and FURIN signals. AlphaGenome predicted the DDX6 regulatory allele downregulates the stress-resistance gene FOXR1 in the nucleus accumbens, while the protective OPRM1 variant (rs1799971) upregulates OPRM1 expression across reward networks. Other signals of interest included IL6R and SHISA9 (EUR); GHR (AFR); and ASTN2 (AMR). Conclusions: This study identifies DDX6 as a novel OUD risk locus, replicates associations with OPRM1 and FURIN, and highlights biologically plausible ancestry-specific signals in AFR and AMR populations. We also replicated top variants in an independent population. Finally, integrating GWAS with deep-learning annotations provides specific, localized biological hypotheses to guide future experimental validation and targeted therapeutics.

Purpose: Nonresponse to routinely collected patient-reported outcome measures (PROMs) threatens the representativeness of aggregated data. We characterized patient-, provider-, and clinic-level factors associated with PROMIS Global-10 nonresponse in routine radiation oncology care. Methods: In this retrospective cohort study, all adults seen at five Mass General Brigham radiation oncology clinics over one year were included. The primary outcome was patient-level nonresponse, defined as never completing the portal-administered Global-10 versus completing it at least once. Using iterative mixed-effects logistic regression, we modeled patient-, provider-, and clinic-level factors. Results: Among 12,214 patients, 71 providers, and five clinics, patient- and appointment-level response rates were 35.4% and 10.9%, with patient-level response ranging nearly fivefold across clinics (12.8% to 66.2%). In Model 1, male sex, lower education, not working, and recent surgery had higher odds of nonresponse, and longer time since diagnosis lower odds. After provider- and clinic-level factors were added, patient sex, education, and employment became nonsignificant, whereas recent surgery (adjusted odds ratio [aOR] 1.97) and longer time since diagnosis (aOR 0.46 for >12 months) persisted. A provider's historical collection rate was protective but attenuated at the clinic level. There, a later program launch (aOR 0.29) and higher historical collection rate (aOR 0.79) correlated with lower nonresponse, whereas academic versus community setting did not. Conclusions: Nonresponse to routinely collected PROMs is a multilevel phenomenon driven substantially by clinic-level implementation factors, not patient characteristics alone. Because response rate is only a proxy for representativeness, PROMs programs and PRO-based performance measures should prioritize representative collection over volume.

Controlled human infection models (CHIMs) are increasingly used in infectious disease research to study pathogen dynamics and evaluate interventions under controlled conditions. However, these studies are resource-intensive and involve ethical and safety constraints, making efficient study design critical. Dose-finding is a key early component in CHIMs, where the aim is to identify a challenge dose that achieves a target infection probability. Traditional rule-based designs are commonly used but can be inefficient, motivating the use of model-based adaptive approaches such as the Bayesian Continual Reassessment Method (CRM). Although CRM has been extensively studied and widely adopted in Phase I oncology trials for identifying the maximum tolerated dose of therapeutics, its application in CHIM settings remains limited, particularly when the endpoint of interest is infection. This tutorial provides step-by-step guidance for implementing a Bayesian CRM in dose-finding CHIMs, using an oropharyngeal Neisseria gonorrhoeae challenge as a motivating case study. The framework outlines key design components, including dose-grid specification, dose-response model, prior elicitation, Bayesian updating, decision rules, and stopping criteria, with particular emphasis on a clinically interpretable parameterisation. Trial operating characteristics are evaluated through simulation studies under multiple dose-response scenarios and prior-predictive analyses, and compared with a commonly used '3+3' type rule-based design. This work highlights the advantages of Bayesian model-based designs for dose-finding in CHIMs over classic rule-based designs and provides a structured, reproducible framework for implementing CRM, supporting their application in future CHIM studies.

Background and Aims The optimal treatment for first episodes and first recurrences of Clostridioides difficile infections (CDI) is unknown and there is emerging evidence for pulse and taper (P-T) regimens. Therefore, we sought to estimate the relative efficacy of treatment options. Methods MEDLINE and CENTRAL were searched from database inception to May 21, 2025 and unpublished conference abstracts were searched from recent infectious disease conferences. RCTs on the treatment of first episodes or first recurrences of CDI comparing fixed-dose or P-T regimens of fidaxomicin or vancomycin were included. The primary and secondary outcomes were 40- and 56-day CDI recurrence, respectively. A random-effects network meta-analysis on the risk ratio (RR) scale was conducted using a standard regimen (10-14 days) of vancomycin as the comparator. Treatments were ranked using the surface under the cumulative ranking curve (SUCRA). Results 8 RCTs were included comprising a total of 2181 patients. For 40-day recurrence, fidaxomicin P-T had the highest probability of ranking best (RR=0.10, 95%Confidence Interval [95%CI]=0.10-0.49, SUCRA=1.00), followed by vancomycin P-T (RR=0.49, 95%CI=0.32-0.76, SUCRA=0.61), fixed-dose fidaxomicin (RR=0.61, 95%CI=0.49-0.76, SUCRA=0.39), and, finally, fixed-dose of vancomycin (SUCRA=0.00). The treatments ranked in the same order for 56-day recurrence, though only 3 RCTs reported on this timepoint. Conclusion Vancomycin P-T, fidaxomicin P-T, and fixed-dose fidaxomicin were all superior to a fixed-dose vancomycin. Head-to-head comparative effectiveness RCTs are needed to quantify their relative effect sizes of and impact on long-term prevention of recurrent CDI.

The emotional content of suicide notes is typically examined using categorical coding, where each labeled passage is treated in isolation from its surrounding language. In contrast, dimensional models of psychopathology propose that affective content varies along continuous gradients. We evaluated this proposition directly. Excerpts from 884 annotated suicide notes were embedded in a semantic space defined solely by their linguistic properties, and we investigated whether human-assigned emotion labels changed smoothly across this space. They did: affective tone showed clear spatial autocorrelation (Moran's $I = 0.18$, $z = 19.68$, $p < 0.001$), an effect that replicated across three different encoders and remained after removing all within-note dependencies. Emotions occupied recognizable yet overlapping regions rather than forming distinct clusters and varied substantially in how tightly they were concentrated: love and hopelessness appeared with similar frequency, but love was far more localized ($z = 15.7$ versus $10.8$). Among all emotions, hopelessness was the most linguistically diffuse, implying that a single categorical label is capturing multiple, qualitatively different manifestations of suicidal distress.

Background The 2022-2023 global mpox outbreak highlighted the critical need for robust genomic surveillance capabilities to track mpox virus (MPXV) evolution and transmission dynamics. Methods Building upon our established SARS-CoV-2 sequencing infrastructure, we implemented a Molecular Loop probe-based long-read sequencing approach using Pacific Biosciences Sequel II technology for comprehensive MPXV genomic surveillance across the United States (US). From August 2024 to June 2025, we generated 326 high-quality whole genome sequences from residual mpox-positive clinical specimens collected by Labcorp across all 10 US Department of Health and Human Services regions. Results Our analysis identified two samples containing clade Ib MPXV in January and June 2025 and captured shifting trends in clade IIb diversity, with 13 distinct lineages observed. We also identified multiple instances of large (~1.6-17.6kb) deletions proximal to the inverted terminal repeats in clade IIb genomes. APOBEC3 mutation analysis indicated substantial evidence of human-to-human transmission among both clades. Further, we observed significantly higher APOBEC3-associated SNPs per kilobase (P<0.001) in clade IIb genomic variable regions relative to their central conserved region. Our assay exhibited strong reproducibility across biological replicates from individual patients and accuracy was confirmed via parallel sequencing of select specimens by US Centers for Disease Control and Prevention (CDC) using metagenomic sequencing. We also demonstrated via custom simulation that our assay discriminates all known MPXV clades and lineages, including those we have not observed in the US. Conclusions Our integrated nationwide surveillance system facilitates real-time genomic tracking of outbreak evolution, with demonstrated capacity across SARS-CoV-2 and MPXV, positioning this platform for rapid deployment during future pathogen emergence.

Background: Monoclonal antibodies have emerged as a prophylactic strategy to prevent symptomatic SARS-CoV-2 infection in immunocompromised individuals. However, the evolutionary and clinical implications of breakthrough infections under this regime remain unclear. Methods: A male in their 80s with a haematological/oncological diagnosis received a 2000 mg intravenous infusion of sotrovimab in March 2023 and was diagnosed with COVID-19 by RT-qPCR from a nasopharyngeal swab in August 2023. Weekly samples (n=24) were collected through February 2024 (171 days). All samples underwent whole-genome sequencing, with select mutations subjected to functional assessment. Findings: Sequencing identified the GE.1 lineage at all timepoints. An intra-host recombination event in ORF1ab (positions 8942-12458) was detected prior to 23 weeks post-detection, followed by a 14-fold increase in viral load (7.42e+06 to 1.00e+08 RNA copies/mL) and a marked shift in the viral population. E340D, a sotrovimab resistance mutation, was detected at low abundance (46%) within the first week post-infection, fluctuated over time, and was nearly fixed by week 15 (107 days) post-detection. We assessed five spike mutations - V36M, S98F, and V213G in the N-terminal domain, Y505P in the receptor-binding domain, and P681Q near the S1/S2 cleavage site - and additionally evaluated the impact of E340D. V36M conferred the highest infectivity across all cell lines, with the most significant effect in low-TMPRSS2 cells. While all mutations showed enhanced infectivity with the addition of E340D, the effect was most pronounced in mutations with lower baseline infectivity. The addition of E340D significantly decreased relative neutralizing titres for V36M, S98F, and V213G, enabling escape from neutralizing antibodies in XBB-responsive individuals, illustrating an enhanced phenotypic advantage. Patient neutralizing activity was absent pre-sotrovimab, and sotrovimab-induced neutralization was further compromised by selection of E340D. Interpretation: Sotrovimab pre-exposure prophylaxis in an immunocompromised patient did not prevent SARS-CoV-2 infection, and selected for resistant mutation E340D, with unexpected fitness consequences across non-receptor binding domain spike regions.

Estimating dengue force of infection (FOI) is essential for understanding transmission dynamics and targeting intervention programmes, yet surveillance data in endemic settings required for estimations are often incomplete, with varying formats. We developed a Bayesian hierarchical catalytic model that jointly fits age-stratified case data, aggregate case data, and seroprevalence surveys within a single framework, incorporating external covariates to improve parameter identifiability. Synthetic validation showed that covariates alone recovered accurate FOI point estimates even when most districts contributed only aggregate data, but did so with poorly calibrated uncertainty; anchoring the model with a single seroprevalence survey was necessary to bring credible interval coverage close to nominal. Applied to 128 districts across Java and Bali, Indonesia (2016-2024), the model revealed substantial spatial heterogeneity in FOI and reporting rates. Many districts in Java exceeded the WHO-suggested seroprevalence threshold for vaccine introduction, yet were classified as low-priority when using reported incidence as prioritisation criterion, particularly in areas with weak surveillance. Model-based seroprevalence estimation, integrating multiple data sources, offers a more consistent basis for identifying high-priority districts for vaccine introduction, and is less susceptible to surveillance bias than reported incidence.

Stochastic epidemic models are a cornerstone of infectious disease epidemiology and are often used to study intervention scenarios. However, large run-to-run variability can make intervention effects difficult to estimate precisely. We revisit the epidemic Sellke construction, which assigns each individual an infection threshold for the cumulative infection hazard such that, conditional on the thresholds, the epidemic trajectory becomes deterministic. This enables coupling of simulations with and without an intervention, yielding low-variance effect estimates even when outcomes such as final size or peak incidence vary widely between runs. We develop an exact, event-driven implementation that maintains infection and recovery events in priority queues. Cumulative infection-hazard updates require O(log N) time per event, yielding overall complexity O(Elog N) for E events in a population of size N. The implementation achieves computational performance comparable to the classical Gillespie algorithm while naturally accommodating non-Markovian infectious periods and complex infectiousness profiles. We illustrate the approach using distance-dependent spread of avian influenza between poultry farms in the Netherlands and a multilayer population with households, schools, and workplaces. In both examples, coupling enables efficient within-run comparisons of intervention scenarios across stochastic realisations.

Patients and caregivers seek informational and emotional support throughout medical care, especially when interpreting unfamiliar laboratory test results. Although resources such as patient portals and online health communities (OHCs) help address questions, gaps remain. The emergence of large language models (LLMs) offers the potential to be a complementary source of support to assist patients and caregivers in understanding and using their test results. The objective of our study is to empirically compare LLM responses to patients online questions containing their laboratory test results to responses written by peers in an OHC. We compared the 519 peer replies to 122 laboratory test-related posts from an OHC to 488 responses generated from four LLMs using mixed computational and qualitative methods. LLMs frequently provided clear explanations of medical terminology and structured interpretations of numeric results but were longer and less readable. Peers offered more personalized, context-specific emotional support. Overall, LLMs have the potential to complement peer responses in OHCs, but require greater emotional depth, reasoning transparency, and alignment with community norms.

Background: Children born prematurely (before 37 weeks) or admitted to the neonatal unit (NNU) are at increased risk of adverse long-term physical health outcomes. It is also recognised that there is an association with later academic performance and special educational needs, however it is not clear whether these broad risk factors could be used as stand-alone heuristics to identify children who may benefit from additional support in educational settings. We aimed to examine the associations between neonatal unit (NNU) admission and educational attainment in mid-childhood. Methods and Findings: Pregnancy data from a prospective birth cohort (Pregnancy Outcome Prediction Study, Cambridge, United Kingdom, 2008-2012) were linked to national educational outcomes (Department for Education, United Kingdom). Multivariable regression models adjusted for maternal, child, and socioeconomic factors were used to evaluate associations between (i) all NNU admissions, (ii) at term NNU admissions >48 hours, (iii) preterm birth without ongoing physical health needs, and educational outcomes at ages 5-11 years. Children who required any NNU care were more likely not to meet expected educational standards across multiple ages and domains in early and mid-childhood: age 5 early year foundation (aOR 1.64, 95% CI 1.19-2.27, p=0.003), phonics at age 6 (aOR 2.43, 95% CI 1.72-3.57, p<0.001), and at age 7 (here assessments were divided into multiple domains): reading (aOR 1.67, 95% CI 1.18-2.38, p=0.004), writing (aOR 1.72, 95% CI 1.25-2.38, p<0.001), mathematics (aOR 1.56, 95% CI 1.09-2.22, p=0.020), and science (aOR 1.85, 95% CI 1.22-2.78, p=0.003). Similar patterns were observed among both at term-born infants who stayed >48hrs in NNU (phonics assessment at age 6 aOR 2.26, 95% CI 1.51-3.36, p<0.001) and in children born preterm without long-term physical health sequelae (phonics assessment at age 6 aOR 3.07, 95% CI 1.96-4.81, p<0.001). These associations were robust to adjustment for demographic, perinatal, and socio-economic factors. By age 11, differences in academic attainment were attenuated and no longer clearly distinguishable across all exposure groups. However, there was an increased likelihood of special educational needs (SEN) at age 11 associated with any NNU admission (aOR 1.78, 95% CI 1.15-2.73, p=0.009), at term NNU admission for >48hrs (aOR 1.88, 95% CI 1.19-3.00, p=0.007), and children born preterm without long-term physical health sequelae (aOR 1.50, 95% CI 1.00-2.25, p=0.049). Predictive performance of any NNU admission for SEN at age 11 was moderate (AUC 0.70, 95% CI: 1.14-2.65, p=0.010), with balanced sensitivity and specificity and high negative predictive value. Conclusions: NNU admission, for both term and preterm infants, is associated with poorer educational outcomes and an increased likelihood of special educational needs in mid-childhood.

Background Post-infection conditions (PICs), such as Long Covid, are associated with heterogeneous, fluctuating symptoms that profoundly affect daily functioning. Despite moderate-certainty evidence from the NIHR-funded LISTEN trial (COV-LT2-0009) that personalised self management support improves outcomes and may reduce societal and economic impacts of Long Covid, many people living with PICs still receive condition-specific services, generic advice, or stand-alone digital tools that do not address their complex needs. Aim To map care approaches in general practice and synthesise UK evidence for PIC management. Design and setting Scoping review and online survey. Method A two-phase study was conducted: (1) a scoping review of UK evidence on PIC management in general practice; and (2) a supplementary online survey of practitioners working in UK general practice to provide contextual insights. Results The scoping review identified 32 studies focused on Long Covid. One study included a comparator group (ME/CFS). Study populations were predominantly white ethnicity and female. Evidence for non-Covid PICs in UK general practice was largely absent. The supplementary survey (n=46) provided preliminary practice-level insights. Healthcare practitioners reported varied PIC presentations, diagnostic uncertainty, limited referral pathways, inequitable access, and low confidence in managing PICs. Conclusion Evidence informing PIC management in UK general practice remains predominantly Long Covid-focused and may not reflect the range of PICs encountered in practice. While survey findings are preliminary and require confirmation in larger samples, they highlight uncertainty around PIC management. Further research is needed to evaluate whether existing Long Covid pathways should be expanded or complemented by broader PIC models. Keywords general practice; Long Covid; self-management; post-viral syndromes